Ravages: An R package for the simulation and analysis of rare variants in multicategory phenotypes.

Current software packages for the analysis and the simulations of rare variants are only available for binary and continuous traits. Ravages provides solutions in a single R package to perform rare variant association tests for multicategory, binary and continuous phenotypes, to simulate datasets under different scenarios and to compute statistical power. Association tests can be run in the whole genome thanks to C++ implementation of most of the functions, using either RAVA-FIRST, a recently developed strategy to filter and analyse genome-wide rare variants, or user-defined candidate regions. Ravages also includes a simulation module that generates genetic data for cases who can be stratified into several subgroups and for controls. Through comparisons with existing programmes, we show that Ravages complements existing tools and will be useful to study the genetic architecture of complex diseases. Ravages is available on the CRAN at https://cran.r-project.org/web/packages/Ravages/ and maintained on Github at https://github.com/genostats/Ravages.

Testing for association with rare variants in the coding and non-coding genome: RAVA-FIRST, a new approach based on CADD deleteriousness score.

Rare variant association tests (RVAT) have been developed to study the contribution of rare variants widely accessible through high-throughput sequencing technologies. RVAT require to aggregate rare variants in testing units and to filter variants to retain only the most likely causal ones. In the exome, genes are natural testing units and variants are usually filtered based on their functional consequences. However, when dealing with whole-genome sequence (WGS) data, both steps are challenging. No natural biological unit is available for aggregating rare variants. Sliding windows procedures have been proposed to circumvent this difficulty, however they are blind to biological information and result in a large number of tests. We propose a new strategy to perform RVAT on WGS data: "RAVA-FIRST" (RAre Variant Association using Functionally-InfoRmed STeps) comprising three steps. (1) New testing units are defined genome-wide based on functionally-adjusted Combined Annotation Dependent Depletion (CADD) scores of variants observed in the gnomAD populations, which are referred to as "CADD regions". (2) A region-dependent filtering of rare variants is applied in each CADD region. (3) A functionally-informed burden test is performed with sub-scores computed for each genomic category within each CADD region. Both on simulations and real data, RAVA-FIRST was found to outperform other WGS-based RVAT. Applied to a WGS dataset of venous thromboembolism patients, we identified an intergenic region on chromosome 18 enriched for rare variants in early-onset patients. This region that was missed by standard sliding windows procedures is included in a TAD region that contains a strong candidate gene. RAVA-FIRST enables new investigations of rare non-coding variants in complex diseases, facilitated by its implementation in the R package Ravages.

Do children with mathematical learning disabilities use the inversion principle to solve three-term arithmetic problems?: The impact of presentation mode.

Numerical inversion is the ability to understand that addition is the opposite of subtraction and vice versa. Three-term arithmetic problems can be solved without calculation using this conceptual shortcut. To verify that this principle is used, inverse problems (a + b - b) can be compared with standard problems (a + b - c). If this principle is used, performance on inverse problems will be higher than performance on standard problems because no calculation is required. To our knowledge, this principle has not been previously studied in children with mathematical learning disabilities (MLD). Our objectives were (a) to study whether 10-year-olds with MLD are able to use this conceptual principle in three-term arithmetic problems and (b) to evaluate the impact of the presentation mode. A total of 64 children with or without MLD solved three-term arithmetic problems (inverse and standard) in two presentation modes (symbolic and picture). The results showed that even though children with MLD have difficulties in performing arithmetic problems, they can do so when the inverse problem is presented with pictures. The picture presentation mode allowed children with MLD to efficiently identify and use the conceptual inversion shortcut and thus to achieve a similar performance to that of typically developing children. These results provide interesting perspectives for the care of children with MLD.

DNA damage in circulating leukocytes measured with the comet assay may predict the risk of death.

The comet assay or single cell gel electrophoresis, is the most common method used to measure strand breaks and a variety of other DNA lesions in human populations. To estimate the risk of overall mortality, mortality by cause, and cancer incidence associated to DNA damage, a cohort of 2,403 healthy individuals (25,978 person-years) screened in 16 laboratories using the comet assay between 1996 and 2016 was followed-up. Kaplan-Meier analysis indicated a worse overall survival in the medium and high tertile of DNA damage (p < 0.001). The effect of DNA damage on survival was modelled according to Cox proportional hazard regression model. The adjusted hazard ratio (HR) was 1.42 (1.06-1.90) for overall mortality, and 1.94 (1.04-3.59) for diseases of the circulatory system in subjects with the highest tertile of DNA damage. The findings of this study provide epidemiological evidence encouraging the implementation of the comet assay in preventive strategies for non-communicable diseases.

Collection and storage of human white blood cells for analysis of DNA damage and repair activity using the comet assay in molecular epidemiology studies.

DNA damage and repair activity are often assessed in blood samples from humans in different types of molecular epidemiology studies. However, it is not always feasible to analyse the samples on the day of collection without any type of storage. For instance, certain studies use repeated sampling of cells from the same subject or samples from different subjects collected at different time-points, and it is desirable to analyse all these samples in the same comet assay experiment. In addition, flawless comet assay analyses on frozen samples open up the possibility of using this technique on biobank material. In this article we discuss the use of cryopreserved peripheral blood mononuclear cells (PBMCs), buffy coat (BC) and whole blood (WB) for analysis of DNA damage and repair using the comet assay. The published literature and the authors' experiences indicate that various types of blood samples can be cryopreserved with only a minor effect on the basal level of DNA damage. There is evidence to suggest that WB and PBMCs can be cryopreserved for several years without much effect on the level of DNA damage. However, care should be taken when cryopreserving WB and BCs. It is possible to use either fresh or frozen samples of blood cells, but results from fresh and frozen cells should not be used in the same dataset. The article outlines detailed protocols for the cryopreservation of PBMCs, BCs and WB samples.

Extension of SKAT to multi-category phenotypes through a geometrical interpretation.

Rare genetic variants are expected to play an important role in disease and several statistical methods have been developed to test for disease association with rare variants, including variance-component tests. These tests however deal only with binary or continuous phenotypes and it is not possible to take advantage of a suspected heterogeneity between subgroups of patients. To address this issue, we extended the popular rare-variant association test SKAT to compare more than two groups of individuals. Simulations under different scenarios were performed that showed gain in power in presence of genetic heterogeneity and minor lack of power in absence of heterogeneity. An application on whole-exome sequencing data from patients with early- or late-onset moyamoya disease also illustrated the advantage of our SKAT extension. Genetic simulations and SKAT extension are implemented in the R package Ravages available on GitHub ( https://github.com/genostats/Ravages ).

Mixed logistic regression in genome-wide association studies.

BACKGROUND: Mixed linear models (MLM) have been widely used to account for population structure in case-control genome-wide association studies, the status being analyzed as a quantitative phenotype. Chen et al. proved in 2016 that this method is inappropriate in some situations and proposed GMMAT, a score test for the mixed logistic regression (MLR). However, this test does not produces an estimation of the variants' effects. We propose two computationally efficient methods to estimate the variants' effects. Their properties and those of other methods (MLM, logistic regression) are evaluated using both simulated and real genomic data from a recent GWAS in two geographically close population in West Africa. RESULTS: We show that, when the disease prevalence differs between population strata, MLM is inappropriate to analyze binary traits. MLR performs the best in all circumstances. The variants' effects are well evaluated by our methods, with a moderate bias when the effect sizes are large. Additionally, we propose a stratified QQ-plot, enhancing the diagnosis of p values inflation or deflation when population strata are not clearly identified in the sample. CONCLUSION: The two proposed methods are implemented in the R package milorGWAS available on the CRAN. Both methods scale up to at least 10,000 individuals. The same computational strategies could be applied to other models (e.g. mixed Cox model for survival analysis).

Adipose Tissue Properties in Tumor-Bearing Breasts.

The tissue stroma plays a major role in tumors' natural history. Most programs for tumor progression are not activated as cell-autonomous processes but under the conditions of cross-talks between tumor and stroma. Adipose tissue is a major component of breast stroma. This study compares adipose tissues in tumor-bearing breasts to those in tumor-free breasts with the intention of defining a signature that could translate into markers of cancer risk. In tumor-bearing breasts, we sampled adipose tissues adjacent to, or distant from the tumor. Parameters studied included: adipocytes size and density, immune cell infiltration, vascularization, secretome and gene expression. Adipose tissues from tumor-bearing breasts, whether adjacent to or distant from the tumor, do not differ from each other by any of these parameters. By contrast, adipose tissues from tumor-bearing breasts have the capacity to secrete twice as much interleukin 8 (IL-8) than those from tumor-free breasts and differentially express a set of 137 genes of which a significant fraction belongs to inflammation, integrin and wnt signaling pathways. These observations show that adipose tissues from tumor-bearing breasts have a distinct physiological status from those from tumor-free breasts. We propose that this constitutive status contributes as a non-cell autonomous process to determine permissiveness for tumor growth.

First genome-wide association study of non-severe malaria in two birth cohorts in Benin.

Recent research efforts to identify genes involved in malaria susceptibility using genome-wide approaches have focused on severe malaria. Here, we present the first GWAS on non-severe malaria designed to identify genetic variants involved in innate immunity or innate resistance mechanisms. Our study was performed on two cohorts of infants from southern Benin (525 and 250 individuals used as discovery and replication cohorts, respectively) closely followed from birth to 18-24 months of age, with an assessment of a space- and time-dependent environmental risk of exposure. Both the recurrence of mild malaria attacks and the recurrence of malaria infections as a whole (symptomatic and asymptomatic) were considered. Post-GWAS functional analyses were performed using positional, eQTL, and chromatin interaction mapping to identify the genes underlying association signals. Our study highlights a role of PTPRT, a tyrosine phosphatase receptor involved in STAT3 pathway, in the protection against both mild malaria attacks and malaria infections (p = 9.70 × 10-8 and p = 1.78 × 10-7, respectively, in the discovery cohort). Strong statistical support was also found for a role of MYLK4 (meta-analysis, p = 5.29 × 10-8 with malaria attacks), and for several other genes, whose biological functions are relevant in malaria infection. Results shows that GWAS on non-severe malaria can successfully identify new candidate genes and inform physiological mechanisms underlying natural protection against malaria.

Rare variant association testing for multicategory phenotype.

Genetic association studies have provided new insights into the genetic variability of human complex traits with a focus mainly on continuous or binary traits. Methods have been proposed to take into account disease heterogeneity between subgroups of patients when studying common variants but none was specifically designed for rare variants. Because rare variants are expected to have stronger effects and to be more heterogeneously distributed among cases than common ones, subgroup analyses might be particularly attractive in this context. To address this issue, we propose an extension of burden tests by using a multinomial regression model, which enables association tests between rare variants and multicategory phenotypes. We evaluated the type I error and the power of two burden tests, CAST and WSS, by simulating data under different scenarios. In the case of genetic heterogeneity between case subgroups, we showed an advantage of multinomial regression over logistic regression, which considers all the cases against the controls. We replicated these results on real data from Moyamoya disease where the burden tests performed better when cases were stratified according to age-of-onset. We implemented the functions for association tests in the R package "Ravages" available on Github.

Detecting the dominance component of heritability in isolated and outbred human populations.

Inconsistencies between published estimates of dominance heritability between studies of human genetic isolates and human outbred populations incite investigation into whether such differences result from particular trait architectures or specific population structures. We analyse simulated datasets, characteristic of genetic isolates and of unrelated individuals, before analysing the isolate of Cilento for various commonly studied traits. We show the strengths of using genetic relationship matrices for variance decomposition over identity-by-descent based methods in a population isolate and that heritability estimates in isolates will avoid the downward biases that may occur in studies of samples of unrelated individuals; irrespective of the simulated distribution of causal variants. Yet, we also show that precise estimates of dominance in isolates are demonstrably problematic in the presence of shared environmental effects and such effects should be accounted for. Nevertheless, we demonstrate how studying isolates can help determine the existence or non-existence of dominance for complex traits, and we find strong indications of non-zero dominance for low-density lipoprotein level in Cilento. Finally, we recommend future study designs to analyse trait variance decomposition from ensemble data across multiple population isolates.

Validation of Gelbond® high-throughput alkaline and Fpg-modified comet assay using a linear mixed model.

Even if the comet assay has been widely used for decades, there is still a need for controlled studies and good mathematical models to assess the variability of the different versions of this assay and in particular to assess potential intra-experimental variability of the high-throughput comet assay. To address this point, we further validate a high-throughput comet assay that uses hydrophilic polyester film (Gelbond®). Experiments were performed using human peripheral blood mononuclear cells (PBMC) either untreated or treated with different concentration of MMS (methyl methanesulfonate). A positive control for the Fpg (Formamidopyrimidine DNA glycosylase)-modified comet assay (Ro 19-8022 with light) was also included. To quantify the sources of variability of the assay, including intradeposit variability, instead of summarizing DNA damage on 50 cells from a deposit by the mean or median of their percentage DNA tail, we analyzed all logit-transformed data with a linear mixed model. The main source of variation in our experimental data is between cells within the same deposit, suggesting genuine variability in the response of the cells rather than variation caused by technical treatment of cell samples. The second source of variation is the inter-experimental variation (day-to-day experiment); the coefficient of this variation for the control was 13.6%. The variation between deposits in the same experiment is negligible. Moreover, there is no systematic bias because of the position of samples on the Gelbond® film nor the position of the films in the electrophoresis tank. This high-throughput comet assay is thus reliable for various applications. Environ. Mol. Mutagen. 59:595-602, 2018. © 2018 Wiley Periodicals, Inc.

Accuracy of heritability estimations in presence of hidden population stratification.

The heritability of a trait is the proportion of its variance explained by genetic factors; it has historically been estimated using familial data. However, new methods have appeared for estimating heritabilities using genomewide data from unrelated individuals. A drawback of this strategy is that population stratification can bias the estimates. Indeed, an environmental factor associated with the phenotype may differ among population subgroups. This factor being associated both with the phenotype and the genetic variation in the population would be a confounder. A common solution consists in adjusting on the first Principal Components (PCs) of the genomic data. We study this procedure on simulated data and on 6000 individuals from the Three-City Study. We analyse the geographical coordinates of the birth cities, which are not genetically determined, but the heritability of which should be overestimated due to population stratification. We also analyse various anthropometric traits. The procedure fails to correct the bias in geographical coordinates heritability estimates. The heritability estimates of the anthropometric traits are affected by the inclusion of the first PC, but not by the following PCs, contrarily to geographical coordinates. We recommend to be cautious with heritability estimates obtained from a large population.

TNFSF10/TRAIL regulates human T4 effector memory lymphocyte radiosensitivity and predicts radiation-induced acute and subacute dermatitis.

Sensitivity of T4 effector-memory (T4EM) lymphocytes to radiation-induced apoptosis shows heritability compatible with a Mendelian mode of transmission. Using gene expression studies and flow cytometry, we show a higher TNF-Related Apoptosis Inducing Ligand (TRAIL/TNFSF10)mRNA level and a higher level of membrane bound TRAIL (mTRAIL) on radiosensitive compared to radioresistant T4EM lymphocytes. Functionally, we show that mTRAIL mediates a pro-apoptotic autocrine signaling after irradiation of T4EM lymphocytes linking mTRAIL expression to T4EM radiosensitivity. Using single marker and multimarker Family-Based Association Testing, we identified 3 SNPs in the TRAIL gene that are significantly associated with T4EM lymphocytes radiosensitivity. Among these 3 SNPs, two are also associated with acute and subacute dermatitis after radiotherapy in breast cancer indicating that T4EM lymphocytes radiosensitivity may be used to predict response to radiotherapy. Altogether, these results show that mTRAIL level regulates the response of T4EM lymphocytes to ionizing radiation and suggest that TRAIL/TNFSF10 genetic variants hold promise as markers of individual radiosensitivity.

Time knowledge acquisition in children aged 6 to 11 years and its relationship with numerical skills.

Acquisition of time knowledge (TK; the correct representation and use of time units) is linked to the development of numerical abilities, but this relationship has not been investigated in children. The current study examined the acquisition of TK and its association with numerical skills. A total of 105 children aged 6 to 11 years were interviewed with our Time Knowledge Questionnaire (TKQ), developed for purposes of this study, and the Zareki-R, a battery for the evaluation of number processing and mental calculation. The TKQ assessed conventional time knowledge (temporal orientation, temporal sequences, relationships between time units, and telling the time on a clock), estimation of longer durations related to birthday and life span, and estimation of the duration of the interview. Time knowledge increased with age, especially from 6 to 8 years, and was strongly linked to numerical skills. Regression analyses showed that four numerical components were implicated in TK: academic knowledge of numbers and number facts (e.g., reading Arabic numerals, mental calculation), number line estimation (e.g., correspondence between a number and a distance), contextual estimation (e.g., many/few leaves on a tree, children in a family), and numerical tasks involving verbal working memory (e.g., comparison of numbers presented orally). Numerical correlations with TK varied according to children's age; subtests based on academic knowledge of numbers, working memory, and number line estimation were linked with TK in the younger children, but only contextual estimation was associated with TK in the older children.

A Bagged, Partially Linear, Tree-Based Regression Procedure for Prediction and Variable Selection.

OBJECTIVES: In genomics, variable selection and prediction accounting for the complex interrelationships between explanatory variables represent major challenges. Tree-based methods are powerful alternatives to classical regression models. We have recently proposed the generalized, partially linear, tree-based regression (GPLTR) procedure that integrates the advantages of generalized linear regression (allowing the incorporation of confounding variables) and of tree-based models. In this work, we use bagging to address a classical concern of tree-based methods: their instability. METHODS: We present a bagged GPLTR procedure and three scores for variable importance. The prediction accuracy and the performance of the scores are assessed by simulation. The use of this procedure is exemplified by the analysis of a lung cancer data set. The aim is to predict the epidermal growth factor receptor (EGFR) mutation based on gene expression measurements, taking into account the ethnicity (confounder variable) and perform variable selection. RESULTS: The procedure performs well in terms of prediction accuracy. The scores differentiate predictive variables from noise variables. Based on a lung adenocarcinoma data set, the procedure achieves good predictive performance for EGFR mutation and selects relevant genes. CONCLUSION: The proposed bagged GPLTR procedure performs well for prediction and variable selection.

Where is the causal variant? On the advantage of the family design over the case-control design in genetic association studies.

Many associated single-nucleotide polymorphisms (SNPs) have been identified by association studies for numerous diseases. However, the association between a SNP and a disease can result from a causal variant in linkage disequilibrium (LD) with the considered SNP. Assuming that the true causal variant is among the genotyped SNPs, other authors demonstrated that the power to discriminate between it and other SNPs in LD is low. Here, we propose to take advantage of the information provided by family data to improve the inference on the causal variant: we exploit the linkage information provided by affected sib pairs to discriminate the causal variant from the associated SNPs. The family-based approach improves discrimination power requiring up to five times less individuals than its case-control equivalent. However, the main advantage of family design is the possibility to carry out the procedure one step further: the linkage information allows inference on causal variants, which are not genotyped but in LD with tag-SNPs displaying association, which is impossible with case-control design. By means of Bayesian methods, we estimate the LD between the observed SNPs and an unobserved causal variant, as well as the allelic odds ratio at the unobserved causal variant. The proposed procedure is illustrated on a multiple sclerosis (MS) family data set including genotypes of SNPs in IL2RA, confirming the advantage of using a family design to identify causal variants. The results of our method on this data suggest the existence of two distinct causal variants in this gene for the MS.

The Use of the Linear Mixed Model in Human Genetics.

We give a short but detailed review of the methods used to deal with linear mixed models (restricted likelihood, AIREML algorithm, best linear unbiased predictors, etc.), with a few original points. Then we describe three common applications of the linear mixed model in contemporary human genetics: association testing (pathways analysis or rare variants association tests), genomic heritability estimates, and correction for population stratification in genome-wide association studies. We also consider the performance of best linear unbiased predictors for prediction in this context, through a simulation study for rare variants in a short genomic region, and through a short theoretical development for genome-wide data. For each of these applications, we discuss the relevance and the impact of modeling genetic effects as random effects.